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We present a universal algorithm for the optimal quantum 
state estimation of an arbitrary finite dimensional system. 
The algorithm specifies a physically realizable positive opera- 
tor valued measurement (POVM) on a finite number of iden- 
tically prepared systems. We illustrate the general formalism 
by applying it to different scenarios of the state estimation 
of independent and identically prepared two-level systems 
(qubits) . 



Suppose we have N quantum objects, each prepared 
in an unknown pure quantum state described by a den- 
sity operator p = |?/;)(?/;|. What kind of measurement 
provides us with the best possible estimation of pi 

Clearly, if we have an unlimited supply of particles in 
state p, i.e. when N approaches infinity, we can estimate 
p with an arbitrary precision. In practice, however, only 
finite and usually small ensembles of identically prepared 
quantum systems are available. This leads to an impor- 
tant problem of the optimal state estimation with fixed 
physical resources. 

A particular version of this problem has recently been 
addressed by Massar and Popescu |l| who have analyzed 
extraction of information from finite ensembles of spin- 
1 /2 particles. They have shown that there exists a funda- 
mental limit on the fidelity of estimation and they have 
constructed a POVM consisting of an infinite continu- 
ous set of operators which attains this limit. It has been 
suggested that the optimal finite POVM may exist, how- 
ever, so far no explicit construction of such finite POVM 
has been provided. We solve this problem by describing 
an algorithm which gives the optimal and finite POVM. 
We also provide a concise derivation of a general bound 
on the fidelity of a state estimation of arbitrary finite- 
dimensional systems (the Massar-Popescu bound follows 
as a particular case). Finite POVMs can be, at least 
in principle, constructed and the optimal state estima- 
tion procedures may be applied in many areas of physics 
ranging from the optimal phase estimation to quantum 
computation. 

In fact, in the following we are solving a more general 
problem of optimal estimation of unitary operations per- 
formed on quantum objects (the state estimation follows 
as a special case). Let us assume that state p is gener- 
ated from a reference state po ~ |?Ao)(V'o| by a unitary 



operation U (x) which is an element of a particular uni- 
tary representation of some group G. Different x denote 
different points of the group (e.g., different angles of ro- 
tation in the case of the SU{2)) and we assume that all 
values of x are equally probable. 

Our task is to design the most general POVM, math- 
ematically described as a set {Or}r=i of positive Hermi- 
tian operators such that '^^r^r — 1 Hi^, which when 
applied to the combined system of all N copies provides 
us with the best possible estimation of p (and therefore 
also of ?7(x)). We quantify the quality of the state esti- 
mation in terms of the mean fidelity 

N times 

r ' " ^ 

/ - V / dx Tr[a C/(x)/5oC/t(x) ®...® J7(x)/5oC/t(x)] 

r JG 

X Tr[i7(x)/5o[/^(x) UrP^Ul], (1) 

which corresponds to a particular choice of a cost func- 
tion used in a context of detection and estimation 
theory. The mean fidelity (|l]) can be understood as fol- 
lows: In order to assess how good a chosen measurement 
is we apply it many times simultaneously on all N parti- 
cles each in state U (x.) poW Ox.) . The parameter x varies 
randomly and isotropically [H over all points of the group 
G during many runs of the measurement. 

For each result r of the measurement, i.e., for each op- 
erator Or, we prescribe the state \ipr) = Ur\tpo) which 
represents our guess (i.e., estimation) of the original 
state. The probability of the outcome r is equal to 
Tr[d^ ;7(x)po;7t(x) ®...® U{x)poW{x)] while the cor- 
responding fidelity of the state estimation is |('0r|''/')P = 
Tr[;7(x)/5o?7t(x) UrPoU}.]. This fidelity is then averaged 
over all possible outcomes and over many independent 
runs of the measurement with randomly and isotropically 
distributed parameters x. We want to find the general- 
ized measurement which maximizes the mean fidelity / 
given by Eq.(0). 

The combined system of N identically prepared refer- 
ence states always remains within the totally symmetric 
subspace of H'' (g) ® . . . , where is fc-dimensional 
Hilbert space of the reference state in which the corre- 
sponding unitary representation J7(x) acts. Thus the di- 
mensionality d of the space in which we construct the 
POVM {Or} isd={ tl^is case the first trace 

in Eq.(|^) can be rewritten as 
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X TV[[/(x)poC/t(x) UrPoU}]d^, (2) 

where J7^(x) is a new representation of the same group 
G; it is equivalent to the A^-fold symmetrized direct prod- 
uct 1^ of the original representation U (x) . Here (x) 
transforms the ( ^j^''~^)-dimensional reference state de- 
noted as r^o- 

We can insert the identity operator U^U^'^ into the 
first trace in Eq. and, taking into account that in 
Eq.(^ we integrate over whole the group G parameter- 
ized by X, we can substitute U^{yi)Uj^^ U^{x) and 
U{x)Uj. U{^)- Now, using the linearity of the trace 
operation as well as the linearity of the representation of 
the group G(UpW is a linear adjoint representation) we 
rewrite Eq.(0) as 

f = Y,Tr[drUr^FUr''% (3) 

r 

where 

F = / C/^(x)f2oC/'^^(x) TV[C/(x)poC/nx) po]rfx, (4) 

JG 

is a positive Hermitian operator. 

Let us now derive an upper bound on the mean fi- 
delity. Taking into account positivity of the operator F 
(i.e., F = J2i^i\'f'i){4'i\i K > 0) and the completeness 
condition for POVM (i.e., J^r = 1) we obtain 

r i r 

< X„,a. U^\C|)^){<|>^\U,''^] (5) 

i r 

r 

From Eq.(^) it clearly follows that the upper bound can 
be achieved if and only if all operators Or forming the 
POVM satisfy the following conditions: 
i) Each Or is proportional to a suitably rotated (by some 
U^) projector on the eigenvector of F with the high- 
est eigenvalue, i.e. for all Or there exists , such that 
Or = cl U^\(j)max){(i}max\Ur^. This , Or morc pre- 
cisely J7r|i/'o)j is our guess associated with the result "r". 
a) All are real and positive, to assure that all O,- are 
positive operators. 

in) Finally, the operators Or have to satisfy the com- 
pleteness criterion Y,r U^\4>max){(iimax\Ur'^ = 1. 

However, these three conditions are not necessarily 
compatible (see Example B below). Therefore, in gen- 
eral, the solution cannot be found via this type of con- 
siderations and we have to proceed in a different way. 



We start with the following observation. Let us as- 
sume that we have some POVM {Or}^^i and the corre- 
sponding guesses Ur which maximize the mean fidelity 
/. We can always construct another POVM with more 
elements which is also optimal. For example, let us con- 
sider a one-parametric subgroup U{(l)) = exjp(iX(j)) of our 
original group G and choose a basis in which 

the action of this subgroup is equivalent to multiplica- 
tion by a factor e*"™''^ (i.e., the operator J7(<^) is diago- 
nal in this basis and ojm are eigenvalues of the generator 
X). Then we take d points (j)a (s = 1, . . . d) and gener- 
ate from each original operator Or a set of d operators 
Ors — ^U'^ {4's)OrU^^<j)s)- In this way we obtain a new 
set of {d ■ R) operators such that the mean fidelity for 
this new set of operators, / = J2r,s T^[Or,s U^sFUj:^s], 
is equal to the mean fidelity of the original POVM {Or} 
because we ascribe to each eventual result [r, s] a new 
guess Ur^s — U{4>s)Ur- Howevcr, in order to guarantee 
that the new set of operators O].^ is indeed a POVM we 
have to satisfy the completeness condition 

1 = E E = E E \u''{^s)drU''\4>s) 

s r s r 

= EE^^E(^0„J->H. (6) 

Let us notice that, by the appropriate choice of (j)g, the 
sum ^-^ — ^ — — can always be made equal to 6m^n 
providing all eigenvalues are non-degenerate |^] (this is 
basically a discrete Fourier transform and we illustrate 
this point in detail in Example A). In this case, the con- 
ditions (||) for the off-diagonal terms in the basis |m) are 
trivially satisfied whereas the diagonal terms are equal 
to unity because the original POVM {O, } guarantees 
that 'Y^^{Or)mm — L Morcover, even if the original set 
of operators {Or} does not satisfy the full completeness 
condition and the conditions for the off-diagonal terms 
are not satisfied (i.e., these operators do not constitute a 
POVM) we can, using our extension ansatz, always con- 
struct a proper POVM {Or.s}- This proves that when we 
maximize the mean fidelity (0) it is enough to assume d 
diagonal conditions rather than the original complete set 
of d^ constraints for diagonal and off-diagonal elements. 

Now we turn back to our original problem of how to 
construct the POVM which maximizes the mean fidelity. 
To do so wc first express the operators Or in the form 
Or = cl U^\^r){^r\U^\ where \^r) are general nor- 
malized states in the d-dimensional space in which the 
operators Or act, and are positive constants. This 
substitution is done without any loss of generality [Q and 
it permits us to rewrite Eq.(^ so that the mean fidelity 
/ does not explicitly depend on C/^, i.e. 

f = Y. clTr[\^r){^r\F]. (7) 
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Obviously, the completeness condition Or 
modified and it reads 



From our discussion above it follows that when maxi- 
mizing the mean fidelity it is enough to apply only d 

constraints | (m|t7Ar rl^'r) | = 1 (here m = 1, . . . , c?) 

out of the (P constraints (^). Therefore to accomplish our 
task we solve a set Langrange equations with d Lagrange 
multipliers Lm- If we express as eigenvalues of the 
operator L — X^m^™ |m)(m| then we obtain the final 
very compact set of equations determining the optimal 
POVM 



From here it follows that \^r) are determined as zero- 
eigenvalue eigenstates. More specifically, they are func- 
tions of d Lagrange multipliers and R vec- 
tors {'^r}r=i [where determine Ur as Ur — [/(x^)]. 
These free parameters are in turn related via R condi- 
tions Det[(F - U^^LU^)] = 0. The mean fidelity now 
is equal to Tri. At this stage we solve a system of d 
linear equations [see the second formula in Eq.(||)] for R 
unknown parameters c^. All solutions for parametri- 
cally depend on L„i and x^ which are specified above. 
We note that the number of free parameters in our prob- 
lem depends on R which has not been specified yet. We 
choose R such that there are enough free parameters so 
that the mean fidelity is maximized and simultaneously 
all are positive. This freedom in the choice of the 
value of R also refiects the fact that there is an infinite 
number of equivalent (i.e., with the same value of the 
mean fidelity) optimal POVMs. The whole algorithm 
is completed by finding (f)g from Eq.(^ which explicitly 
determine the finite optimal POVM {0'^^]. This is the 
main result of our Letter. 

In the following we will apply this general algorithm 
into two simple examples: 
Example A 

Suppose we have N identical copies of spin 1/2 all pre- 
pared in the same but unknown pure quantum state. 
If we chose the group G to be C/(2), i.e. the complete 
unitary group transforming a two-level quantum system, 
we can straightforwardly apply the optimal estimation 
scheme as described above. To be more precise, due to 
the fact that there exist elements of the group U{2) for 
which the reference state is the fixed point (i.e., it is in- 
sensitive to its action) we have to work only with the 
coset space ^^^^^u(n~i) Hi- In the present case this is 
a subset of the SU{2) group parameterized by two Euler 
angles 0, ijj (the third Euler angle x is fixed and equal to 
zero). This subset is isomorphic to the Poincare sphere. 

The unitary representation U is now the representa- 
tion (i) (we use a standard classification of SU{2) rep- 
resentations, where (j) is the spin number). Its iV-fold 



1 is now symmetrized direct product (we denote this representa- 
tion as U^) is the representation classified as (^) (which 
transforms a spin-A''/2 particle). Choosing the standard 
basis \j, m) with m = —j, . . . j in which the coordinate 
expression for U{9,ip) corresponds to standard rotation 
matrices DJ„_„(6l, V', 0) = c-""'^ d^,„(6') §, we obtain 
the matrix expression for the operator F 



(8) 



p — 



27r 



sm{d)dd 



8tt 



(1 



COS ( 



(10) 



X D' 



D 



N/2 



{N + 2){N + 1) 



When we insert this operator in the Eq.(^ we immedi- 
ately find the upper bound on the mean fidelity to be 
equal to This is the result derived by Massar and 

Popescu |1| who have also shown that this upper bound 
can be attained using the special POVM which consists 
of an infinite continuous set of operators proportional to 
isotropically rotated projector j^, t)(T' tI S- 

However, following our algorithm, we can now con- 
struct the optimal POVM which is finite. To do so, we 
have to find a finite set of pairs of angles {{6r,4'r)} such 
that the completeness conditions (^) which now take the 
form 

J2cl e-''^'-('"-")# „(0.) dl^iOr) = <5™,„, (11) 



are fulfilled. Following our general scheme we first sat- 
isfy the completeness conditions ( p]) for diagonal terms 
[compare with Eq.(0)] 



= 1; 



= -N/2, ...N/2. (12) 



To satisfy these completeness conditions we choose + 1 
angles 9r to be equidistantly distributed in the (0, tt) (ob- 
viously, there are many other choices which may suite the 
purpose ~ see discussion below Eq.(||)). Then we solve 
the system of linear equations for A^ -|- 1 variables . For 
this choice of 9r the system (^2|) has non-negative solu- 
tions. Finally we satisfy the off-diagonal conditions by 
choosing A^ -I- 1 angles V-'s = for each 9r. In this 

^i^^y ^ ^^^^ y ^ -N/2,... N/2 



N+ 



+ 1 Z^s=0 ' 



and the off-diagonal conditions are satisfied straightfor- 
wardly. This concludes the construction of the optimal 
and finite POVM for the spin-1/2 state estimation. 
Example B 

The algorithm can also be used to estimate a unitary 
evolution of quantum systems. Consider, for example, 
a system of A^ qubits, all initially prepared in some ref- 
erence state and undergoing a free evolution described 
by the U[l) group (e.g. like in the Ramsey type exper- 
iments pOl). Our task is to find a measurement which 
provides the optimal estimation of the phase of the U{1) 
rotation. N.B. the phase estimation without any a pri- 
ori information is different from the frequency standards 
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experiments, where the issue is the abihty to distigu- 
ish neighbouring quantum states with the best resolution 
(see, for example pl). 

All unitary irreducible representations of U{1) are one 
dimensional and they are parameterized by a single inte- 
ger number (j). Rotation of a single isolated qubit is de- 
scribed by a representation U classified as (0) ® (1) (here 
© denotes a direct sum of representations). The entire 
system of TV qubits (these are assumed to be nonentan- 
gled) is then transformed due to the representation 
which is equal to the symmetrized A^-fold direct product 
of the basic representation U. The representation is 
equal to a direct sum of irreducible representations of the 
form (0)® (l)e. . . {N) which act in the -1-1 dimensional 
space spanned by basis vectors \m), m = 0, 1, ... TV. In 
this basis matrix elements F^.n of the operator F given 
by Eq.(|[) take the form 



Ffn.n — 



^'^ d^J \J (]v-m)(]v-n) 



27r 



2N+1 



g,:V.(n-m) (1+cOS^) 



\J (jV-m)(jV-n) 



). (13) 



The upper bound on the fidelity Eq. (|^) is now too con- 
servative to be of any use (greater than unity). We can, 
however, solve the system of Eqs. (^ which in this par- 
ticular case reads 



F -L I*) = 0, 



|(m|«')p = l; Vm. 



The condition Det(i^ — L) — now determines the eigen- 
vector 1 5') with the zero eigenvalue as a function of La- 
grange multipliers L„j. When we substitute this eigen- 
vector into the second equation in Eq.(|lJ) we obtain a 
set of equations for Lm from which the reference state 
1 4') can be determined. The final POVM is then con- 
structed by rotation of 1 5*) by iV -|- 1 angles 0s in such a 
way that all off-diagonal elements oi^^{Os)m,n become 
equal to zero. This is done in exactly the same way as 
in Example A. The resulting POVM corresponds to the 
von Neumann measurement performed on the composite 
system of all N ions. 

We note that the mean fidelity / has a rapidly growing 
number of local extrcma which originate from the roots of 
certain polynomials. From this set of extrema we easily 
choose the global maximum which corresponds to the 
optimal POVM. 

The maximal mean fidelity for the first six N 

3 2+V2 ll+2\/3 
4' 4 ' 16 ' 



is: 



32+10v^-l-3v^+\/6 
64 



5+\/6 32+5\/2+2\/5+5V2\/ 3+2x/2 
8 ' 64 ' 

. . which approximately gives: 0.750, 



0.854, 0.904, 0.931, 0.947, 0.957, . . . 

In conclusion, we have presented a general algorithm 
for the optimal state estimation from finite ensembles. 
It provides finite POVMs which, following the Neumark 
theorem p^ , can, at least in principle, be implemented 
as simple quantum computations. 
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